Cluster Dynamical Mean Field Theory of the Mott Transition 
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We address the nature of the Mott transition in the Hubbard model at half-filling using cluster 
Dynamical Mean Field Theory (DMFT) . We compare cluster DMFT results with those of single site 
DMFT. We show that inclusion of the short range correlations on top of the on-site correlations, 
already treated exactly in single site DMFT, do not change the order of the transition between the 
paramagnetic metal and the paramagnetic Mott insulator, which remains first order. However, the 
short range correlations reduce substantially the critical U and modify the shape of the transition 
lines. Moreover, they lead to very different physical properties of the metallic and insulating phases 
near the transition, in particular in the region of the phase diagram where the two solutions coexist. 
Approaching the transition from the metallic side, we find an anomalous metallic state with very 
low coherence scale at temperatures as low as T — O.Olt. The insulating state is characterized by 
the relatively narrow Mott gap with pronounced peaks at the gap edge. 



The correlation driven metal insulator transition is one 
of the most fundamental problems in condensed matter 
physics, and continues to receive intensive attention. It 
is realized in numerous transition metal oxides and some 
organic salts, by application of the pressure or isovalent 
chemical substitutions [![. The metallic state far from 
the transition is well described by the Fermi liquid the- 
ory, illustrating the wave-like properties of electrons in 
solids. In the insulating side, the electron behaves as 
a localized particle. Near the transition, the effective 
Coulomb repulsion between the carriers is of the same 
order as the kinetic energy term in the Hamiltonian. 
This regime probes the dual character of electron, namely 
the particle- and wave-like character, and requires a non- 
perturbative method for its description. 

The nature of the metal to insulator transition depends 
strongly on the degree of magnetic frustration. In the 
limit of very large magnetic frustration, the insulating 
state is a simple paramagnetic state with local moments 
carrying log(2) entropy. The metallic state is a Fermi liq- 
uid with a very heavy mass. The mass increases as the 
transition is approached to match the large entropy of the 
frustrated paramagnetic insulator. This is the essence of 
the Brinkman-Ricc theory of the metal insulator transi- 
tion, which has been substantially extended by the single 
site DMFT of the Hubbard model in the paramagnetic 
phase • The key predictions of this approach, such as 
the existence of a first order line ending in a second order 
Ising point, and numerous high temperature crossovers, 
have been verified experimentally [3(. The first order 
phase transition in a strongly frustrated situation has 
been confirmed by cluster DMFT studies [H, 0] and by 
other techniques [a]. 

The completely unfrustrated case is also well under- 
stood along the lines first drawn by Slater, and realized 
in the half filled one band Hubbard model with only near- 
est neighbor hoppings. Here, the metal insulator transi- 
tion is driven by the long range magnetic ordering. The 
system is insulating and magnetic for arbitrarily small 



values of U, as a reflection of the perfect nesting of the 
band structure. The insulating gap results from the for- 
mation of a spin density wave that Bragg scatters the 
electronic quasiparticles. 

The character of the metal insulator transition with an 
intermediate degree of frustration (when the long range 
magnetic order is fully suppressed, but with strong short 
range magnetic correlations) remains an open problem. 
Qualitative modifications of the character of the transi- 
tion are expected, since at low temperatures the para- 
magnetic insulating state has very low entropy. This 
problem can be addressed by a sharp mathematical for- 
mulation studying the paramagnetic solution of the clus- 
ter DMFT equations of the Hubbard model, keeping the 
short range correlations only. The early cluster DMFT 
studies received conflicting answers depending on differ- 
ent cluster schemes and different impurity solvers 0, HI • 
However, by going to very low temperatures using new 
algorithmic developments, we completely settle this ques- 
tion. 

Method: To study this problem, we apply cellular dy- 
namical mean field theory (CDMFT) [3 to the two- 
dimensional Hubbard model, using plaquette as a ref- 
erence frame. In this formalism, the lattice problem is 
divided into 2x2 plaquettes and the lattice problem is 
mapped to an auxiliary cluster quantum impurity prob- 
lem embedded in a self-consistent electronic bath. The 
latter is represented by an 8 x 8 matrix of impurity hy- 
bridization A, which is determined by the condition 



A(iu>) = icu + fx — Y> c (iLo) — 



E- 



i 



iuj + it, — t c (k) — Y> c (iui) 



where S c is the matrix of cluster self-energies, t c {k) is 
the matrix of tight-binding hoppings expressed in terms 
of the large unit cell (2 x 2) of the cluster, and k runs over 
the reduced Brillouin zone of the problem. We choose 
the two dimensional square lattice with only the nearest 
neighbor hopping t. 
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The cellular DMFT approach has already given nu- 
merous insights into frustrated models of kappa organ- 
ics 0, El as well as the doping driven Mott transition 
in the Hubbard model, when treated with a variety of 
impurity solvers 12j. In this letter, the auxiliary cluster 
problem is solved with the numerically exact continuous 
time quantum Monte Carlo (CTQMC) method (I3.[l4|. 
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FIG. 1: (Color online) a) The phase diagram of the param- 
agnetic half-filled Hubbard model within plaquette-CDMFT. 
Inset: The histogram of the two insulating states. It shows 
the probability for a given cluster eigenstate among the 16 
eigenstates of the half-filled plaquette. The singlet plaquette 
ground state has the highest probability, b) For comparison, 
the corresponding phase diagram of the single site DMFT 
(using the same 2D density of states) is shown. The coexis- 
tence region is shown as the shaded region. The dashed line 
marks the crossover above the critical point. The crossover 
line was determined by the condition that the imaginary part 
of the self-energy at few lowest Matsubara frequencies is flat 
at the crossover value of U. For easier comparison, the x-axis 
is rescaled and the reduced value of U r = u ~ Umit is used. 

U MIT 

The critical value of U is Umit = 6.054 in cluster case and 
Umit = 9.35t in single site case. Pentagons in panel a) mark 
the points in phase diagram for which we present the local 
spectral functions in Fig. 2. 



Results: FigOJi shows the phase diagram of the Hub- 
bard model within cluster DMFT at half-filling in the 
absence of long range order. For interaction strength 
U < U C 2(T), we find a metallic solution while for U > 
U c i(T), a Mott insulating solution exists. The two tran- 
sition lines Ud(T) and U C 2{T) cross at a second order 
cndpoint, at temperature Tmit ~ 0.09t and interaction 
strength U mit ~ 6.054. It is clear that one of the most 



salient features of the single site DMFT phase diagram 
(shown in Fig. [Tp), namely the existence of a first order 
phase transition, survives in plaquette-DMFT. 

Still there are substantial modifications to the sin- 
gle site DMFT results when U/t is close to its critical 
value. Namely, i) Strong short ranged antiferromagnctic 
correlations significantly reduce the value of critical U 
at which the second order endpoint occurs. Note that 
the plaquette-DMFT critical C/(~ 6.05i) is in very fa- 
vorable agreement with the Monte-Carlo crossover U at 
which the pseudogap develops at intermediate temper- 
atures accessible by determinantal Monte Carlo (figure 



Ref. 



15l ). This critical U will increase if the sys- 
tem is more frustrated at short distance. For example, 
the inclusion of the next nearest hopping t' has this ef- 
fect and was studied in Ref. [lj| ii) The shape of the 
coexistence region, where both metallic and insulating 
solutions exist, is significantly different. The high tem- 
perature crossover lines (dashed line above T ~ O.lt in 
FigfT]) are similar since at high temperature the entropy 
of the paramagnetic insulator is of the order of log(2) in 
both cluster and single site approach. As the tempera- 
ture is increased, the large entropy insulating state wins 
over the lower entropy metallic state. At low tempera- 
ture, the situation is very different. In single site DMFT, 
the metal wins at low temperature in the transition re- 
gion because the emergence of the itinerant quasiparticle 
inside the Mott gap lowers the free energy of the strongly 
disordered Mott state. In the cluster case, the Mott in- 
sulator at very low temperature is very different and has 
small entropy due to short range singlet formation. The 
small entropy of this state can be confirmed by the "va- 
lence histogram" shown in the inset of Fig. la. The high 
temperature insulating state, which has entropy of the 
order of log(2), populates many states of the plaquette 
with significant probability. In contrast, there is only one 
significant eigenvalue of the density matrix in low tem- 
perature, corresponding to the singlet state. The insulat- 
ing phase at low temperature has thus very small entropy, 
and the bad metal has larger entropy, hence decreasing 
temperature favors insulator over metal. The actual first 
order line (dashed line in Fig. [IJi inside the coexistence 
region, where the free energy of the two phases equals) 
therefore bends back and critical U decreases with de- 
creasing temperature. It is apparent that the zero tem- 
perature transition in cluster-DMFT happens at U c \ and 
not at U c2 as in DMFT. 

While the shape of the DMFT phase diagram strongly 
resembles the phase diagram of the Cr-doped V2O3, the 
reentrant shape of the cluster-DMFT transition resem- 
bles more the k— organic diagram 17] as pointed out in 
Ref. |. 

To understand the effects brought about by the short 
range magnetic correlations near the transition, we fo- 
cus on the local spectral functions displayed in Fig. [5] 
As in single site DMFT, below U& (Fig. [3^) the system 
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FIG. 2: (Color online) The local spectral function for four rep- 
resentative values of U/ts and temperature T = O.Olt marked 
by pentagons in Fig.l. a) For U below U c i the system is 
in Fermi liquid regime with rather large coherence tempera- 
ture, b) In the coexistence region, the insulating solution has 
a small but finite gap (~ 0.2t). c) The metallic solution in 
the same region is strongly incoherent and the value at zero 
frequency decreases due to the finite scattering rate (see self- 
energy in Fig.|3K). d) For U above U C 2, the Mott gap steadily 
increases with U. 
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is a normal Fermi liquid with a reduced width of the 
quasiparticle peak (Z ~ 0.4) and well developed Hub- 
bard bands around —2.5t and 2.5i. 

The insulator in the coexistence region (Fig. [Jb) is 
however very different than Mott insulator in single site 
DMFT. The Mott gap is small and it vanishes at U c i 
where the insulating solution ceases to exist. At low tem- 
perature very pronounced peaks at the gap edge appear. 
These peaks are a clear hallmark of the coherence peaks 
characteristic of a Slater spin density wave. This has 
been noticed earlier in numerous studies of the Hubbard 
model d[H, [3, as well as in the single site DMFT so- 
lution in the ordered phase of the unfrustratcd lattice, 
which captures the physics of perfect nesting. 

With increasing U above U C 2 (Fig. [2H) the Mott gap 
increases but the peaks at the gap edge remain very pro- 
nounced. Only at very large U comparable to the crit- 
ical U of the single site DMFT they lose some of their 
strength and dissolve into a featureless Hubbard band. 

The metallic state, which competes with the insulator 
in the coexistence region, (Fig. has similar width of 
the quasiparticle peak as the Fermi liquid state at U < 
U c i. Hence the quasiparticle renormalization amplitude, 
as extracted at finite but low temperature T = O.Olt is 
rather large. On the other hand, this metallic solution 
has somewhat reduced height of the quasiparticle peak 
which is mostly due to incoherent nature of the solution. 

The incoherence can also be identified from the raw 



FIG. 3: (Color online) top: The imaginary part of the clus- 
ter self energies for the same parameters as in Fig. 2. Due 
to particle-hole symmetry, the (n, n) and (0, 0) cluster self- 
energies have the same imaginary part and we show only one 
of them. Below the metal-insulator transition shown here in 
a), the momentum dependence of the self-energy is rather 
weak and the cluster solution is very similar to the single site 
DMFT solution. In the coexistence region, the metallic so- 
lution shown here in c) is strongly incoherent especially in 
the (it, 0) orbital. For the insulating solutions in b) and d), 
the (tv, 0) scattering rate diverges which opens the gap in the 
spectra, bottom: e) _RcEk(0 _ ) — /j, as a function of U. Due 
to particle-hole symmetry, i?eE K= (^ j0 ) — A* vanishes. 



data on the imaginary axis. In Fig. [3] we show the imag- 
inary self-energy for the different cluster momenta if, 
which can be thought as the orbitals of the multi-orbital 
model associated with the cluster. In plaquette geometry, 
the self-energy is diagonal in cluster momentum base and 
the on-site, nearest-neighbor, and next-nearest-neighbor 
self-energies can be constructed as the linear combination 
of these orbital self-energies []~2| . 

Below U c i , the self-energies of all four orbitals are very 
similar and results are close to the single-site DMFT. 
The metallic phase in the coexistence region Fig. [3}; has 
a large scattering rate in the (tt, 0) orbital, in the orbital 
which contributes most of the spectral weight at the fermi 
level. The coherence scale in this strongly incoherent 
metal is thus severely reduced. The scattering rate as 
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a function of temperature is not quadratic even at T = 
0.0 li and remains large ~ 0.2t at that temperature. 

In Fig. [Sb, d the Mott insulating state can be identi- 
fied by the diverging imaginary part of the 
Due to particle hole symmetry, the real part of the same 
quantity vanishes. Therefore, the only way to open a 
gap in the single particle spectrum is to develop a pole 
at zero frequency £(^,0) — C/(iw). We checked that the 
insulating state in the coexistence region has the char- 
acteristic l/{iuj) behavior at very low temperature and 
the coefficient C in the coexistence region decreases as U 
decreases. The closure of the gap at the U c \ transition 
point is confirmed by the vanishing of C at that point. 

The other two orbitals expel their Fermi surfaces by a 
different mechanism identified in Ref . [2(3, namely the real 
parts of the self-energy are such that the effective chem- 
ical potential jU e // = A* — ^(0~) moves out of the band. 
The separation of the two orbitals gradually increases as 
U increases, and it jumps at the critical U showing the 
hysteresis behavior displayed in Fig. [3^. 




iif=(7r,0),(0,Jr) 
»if=(0,0),(7r,jr) 
DMFT 



We therefore can not determine the low energy Z which 
might vanish at U C 2 at zero temperature. 

In conclusion, we used essentially exact numerical 
method, continuous time quantum Monte Carlo, and 
clarified the nature of the Mott transition in plaqucttc- 
DMFT. The short range correlations which are accounted 
for in this study but are absent in single site DMFT do 
not change the order of the Mott transition which re- 
mains first order with coexistence of metallic and insu- 
lating solution. Our cluster DMFT study predicts the 
existence of anomalous metallic state within the coexis- 
tence region with very low coherence temperature. This 
regime could be relevant to the interpretation of experi- 
ments in VO2 2l| and PrNi03 under the applied pres- 
sure j22| where an anomalous metallic state was reported. 
On the theoretical side, the plaquette DMFT brings new 
light on the nature of the interaction driven MIT. The 
cluster DMFT of this problem retains aspects of Mott 
physics, as described in single site DMFT, and Slater 
physics. It does that by having two orbitals ((n, 0) and 
(0, 7r)) exhibit a Mott transition while the remaining or- 
bitals ((0,0) and (71", 7r)) undergo a band transition. This 
Slater-Mott transition requires momentum space differ- 
entiation and has no analog in single site DMFT. 
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FIG. 4: (Color online) The quasi-particle residue Z vs U/t 
for different orbitals in CDMFT. Below the transition point, 
the (0, 0) and (tt, tt) orbitals have essentially the same Z as 
single site DMFT (dotted line) while the quasiparticles are 
more renormalized in (n, 0) orbital. 

The important issue in the metal insulator transition 
(MIT) is whether the short range magnetic exchange in 
the Hubbard type of models allow the Brinkman-Rice 
scenario of diverging effective mass. In Fig. U] we plot the 
quasiparticle renormalization amplitude Z of the four dif- 
ferent orbitals of the plaquette. As shown in Fig. QJ the 
growth of the effective mass in cluster DMFT is cut-off 
by the exchange interaction and the spatial coherence is 
lost way before the quasiparticles acquire a large effective 
mass. The lattice Z^ is a linear combination of the two 
values plotted in Fig. 0J The quasiparticles at (tt, 0) and 
(0, 7r) are renormalized more strongly than those away 
from the two points. More importantly, close to U&, 
where the system is still coherent at T = i/100, the quasi- 
particle renormalization amplitude is rather large for the 
plaquette without frustration [Z ~ 0.36). Very near and 
inside the coexistence region, the metallic state remains 
very incoherent at our lowest temperature T = i/100. 
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